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Abstract 

We explicitly construct and study the statistics of flux vacua for type IIB 
string theory on an orientifold of the Calabi-Yau hypersurface Ft 1 2 2 6 j , 
parametrised by two relevant complex structure moduli. We solve for 
these moduli and the dilaton field in terms of the set of integers defin- 
ing the 3-form fluxes and examine the distribution of vacua. We compare 
our numerical results with the predictions of the Ashok-Douglas density 
det(— 1Z — u), finding good overall agreement in different regions of moduli 
space. The number of vacua are found to scale with the distance in flux 
space. Vacua cluster in the region close to the conifold singularity. Large 
supersymmetry breaking is more generic but supersymmetric and hierar- 
chical supersymmetry breaking vacua can also be obtained. In particular, 
the small superpotentials and large dilaton VEVs needed to obtain de Sit- 
ter space in a controllable approximation are possible but not generic. We 
argue that in a general flux compactification, the rank of the gauge group 
coming from D3 branes could be statistically preferred to be very small. 
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1 Introduction 

It has been claimed that the current status of string theory resembles that of par- 
ticle physics in the early 1960's[l j in the sense that there are many experimental 
results but no organising principle. For string theory the place of the experiments 
is taken by the many explicit string vacua that have been constructed over the 
years. One desirable avenue towards understanding string theory involves the 
development and questioning of its basic principles. Another less ambitious but 
equally important program is to explore the large variety of string vacua with 
the hope of learning more about the structure of the theory and of finding some 
of its phenomeno logical and/or cosmo logical implications. 

In the phenomenological direction there has been substantial progress through 
the construction of string vacua with properties very similar to the standard 
model^. In cosmological applications, concrete ways have been found to realise 
inflation within string theory ^j. Furthermore the construction and study of 
string vacua and their properties has played an important role over the years 
in the understanding of the theory itself, as exemplified by the discovery of T- 
duality [4j, mirror symmetry [^j and S'-duality 0. Mirror symmetry provides 
probably the most direct example in the sense that even though there were sev- 
eral ideas about its existence, the actual realisation came only after a plot of 
many of the vacua was found to be essentially symmetric under the exchange of 
complex and Kahler structure moduli, thus providing 'experimental' evidence for 
this symmetry. 

Whether we like it or not, the lack of a close experimental test of the theory 
and the large abundance of vacua suggest that a reasonable direction to pur- 
sue is to create a large database of string vacua as a temporary substitute for 
experimental results. This database may eventually be used to extract model- 
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independent features of string vacua that may be testable. It may also give us 
some further hints into the structure of the theory itself. 

These issues have been recently clarified through the study of antisymmetric 
tensor flux compactifications IB]. In the absence of these fluxes there were 
myriads of vacua, each with a large continuous degeneracy parametrised by the 
moduli of the compactification space. Turning on the fluxes substantially reduces 
this degeneracy by making it discrete, owing to the Dirac quantisation condition 
for the fluxes. The fluxes can fix the values of the moduli, which then belong to 
a discrete set determined by the allowed fluxes. This huge discrete degeneracy 
makes a statistical analysis suitable for the study of vacua. The stabilisation 
of the dilaton and complex structure moduli is relatively straightforward, but 
flux compactifications may also enable the construction of metastable vacua with 
positive cosmological constant with the Kahler structure moduli also stabilised 
as in the KKLT scenario [01 E]] . 

Even though the number of flux vacua appears to be enormous for a general 
Calabi-Yau manifold, there are a number of requirements that must be satisfied 
for these vacua to be acceptable. First, as the entire framework is founded on a 
weak coupling approximation, only solutions leading to values of Im r ^> 1, where 
r = Co + ie~^ is the dilaton-axion field that measures the string coupling. The 
self-duality of type IIB under the SL(2, Z) ^-duality symmetry implies that any 
solution can be mapped to the fundamental domain. Having mapped solutions 
to this region we can then only trust values of the dilaton substantially greater 
than one. As we will discuss in section 12.21 this is expected to be naturally 
suppressed statistically, with the proportion of vacua with Imr > tq falling as 
j^^. Secondly, the fluxes may or may not break supersymmetry by themselves 
and it is of interest to know the value of the supersymmetry breaking scale for 
each model. This is more urgent if the models are embedded in the full KKLT 
scenario in which non-perturbative superpotentials are also included to fix the 
Kahler moduli 3 . In order to guarantee a volume sufficiently larger than the string 
scale and believe the four- dimensional effective field theory analysis, the value of 
the gravitino mass induced only by the fluxes has to be hierarchically small, a 
condition that is not straightforward to satisfy. 

Recently, Douglas and collaborators have developed statistical techniques to 
study string vacua e.g. jT^l EH QZj- I n particular, in [Tfij . a general formula was 
proposed that captures the statistical nature of the vacua, in terms of a density 
given by det(— 1Z — to) where TZ is the curvature two form and u the Kahler two- 
form on moduli space. In [TH] the validity of this formula was successfully tested 
in a one-modulus example. 

3 In the KKLT scenario, the combination of fluxes that themselves break supersymmetry 
with a non-perturbative superpotential results in a supersymmetric minimum. In this case 
supersymmetry is broken only after the introduction of anti D3-branes [51 111), magnetic fluxes 
on D7 branes by a' corrections or m other vacua of the complex structure moduli |T3|. 
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In this note we extend this analysis to a simple two-modulus compactifica- 
tion, corresponding to the three-fold realisable as a hypersurface in the weighted 
projective space P[i,i,2,2,6]- We turn on 3-form fluxes and use the techniques of 
Candelas et al (THl 1201 I2J] to compute the periods. This allows us to determine 
the superpotential for the dilaton and complex structure moduli fields. We use 
a Monte Carlo analysis to study the distribution of vacua in the vicinities of the 
Landau-Ginzburg and conifold regions of moduli space. We compute the Ashok- 
Douglas index and compare with our results, finding good general agreement. 
We also investigate the distribution of values of the dilaton field as well as the 
supersymmetry breaking scale. 

While in the last stages of this project we received an article |22j that dis- 
cusses the Ashok-Douglas density for the same Calabi-Yau without performing 
the numerical analysis. 

2 General Framework 

We here review the construction of the particular Calabi-Yau we study, and the 
calculation of its periods[2IIJ Some explicit flux compactifications on this 
Calabi-Yau are also presented in [25] . 

2.1 The Model 

There exist many Calabi-Yau manifolds that can be realised as hypersurfaces in 
weighted projective space. The weighted projective space PjL ki k2 ki fc i is defined 
by 

Oo, w 1} w 2 , w 3 , w 4 ) = (X k °w , X kl wi, X k2 w 2 , X k3 w 3 , X ki w 4 ) 

and has four complex dimensions. To obtain a space of three complex dimensions 
we restrict to the hypersurface P{wj) = 0, where P is a polynomial in the u^. 
The condition that such a hypersurface be Calabi-Yau is that deg(P) = Ylt=o ki- 
ln this paper we will make use of the Calabi-Yau hypersurface in Ft x 2 2 6 j, which 
arises from 

wl 2 + w 12 + w% + wt + wt = 0. (1) 

We shall denote this manifold by M. M has h 1 ' 1 = 2 and h 2 ' 1 = 128. The 
number of complex structure moduli is determined by the number of monomial 
deformations of degree 12 that can be added to (0). There are two complex 
structure moduli that are particularly important, which we shall denote by tp 
and 0. These perturb ((TJ) to 

P(v)i) = wl 2 + wl 2 + wl + + w\ - 12ip - 20 (wowf) =0. (2) 

Their significance is due to their survival in the mirror manifold W, obtained by 
identifying points in M. related by the action of G, the maximal group of scaling 
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symmetries of (JTJ. G here is Z 2 x Zg, and its action is represented by 

Z 2 : (wi,w 4 ) — »■ (a^i, awi) with a 2 = 1, 

Z 6 : (w 1 ,w 3 ) -> (/3V,« with /3 6 = 1, 
Z 6 : (tUi,u> 2 ) — * (7 5 Wi,7iy 2 ) with 7 s = 1. 

These manifestly leave equation (j2j) invariant. All other degree 12 deformations 
of equation (J2J) are not well defined on the mirror. As W has ft, 1 ' 1 = 128 and 
ft, 2,1 = 2, it is an example of a Calabi-Yau with two complex structure moduli. 
Both M. and W develop a conifold singularity when 864^ 6 +0 = 1. This condition 
follows from requiring P(wi) = VP(u>i) = with V VP non-singular. 

The Ashok-Douglas index density describes the statistics of the class of mod- 
els described in [E], namely flux compactifications of type IIB orientifolds. To 
fit the above model into this framework, we embed it into an F-theory model 
with an orientifold limit according to Sen's prescription [24J. The correspond- 
ing elliptically fibered four- fold is the hyper surf ace in weighted projective space 
Pm 1 2 2 12 18 j with x — 19728. The associated tadpole cancellation condition then 
reads [23]: 

A^ds + Afflux = y a = 822 = L - ( 3 ) 

Here is the net number of D3 branes and iVfl ux the contribution of the fluxes 
to the D3 brane charge, 

^ = 7^T 2 I H ?J A F 3 . (4) 
(2tc) a' 2 Jm 

Here H s = dB 2 and F 3 = dC 2 are the NS-NS and R-R 3-forms respectively. 
To cancel tadpoles in a supersymmetric fashion we then require the fluxes to 
carry at most 822 units of D3-brane charge (this number as usual measures the 
amount of negative D3 brane charge coming from D7 branes and orientifolds, in 
the orientifold realisation of the model). 

In the low-energy approximation, the moduli of the Calabi-Yau appear as 
scalar fields in a supergravity theory. The effective four- dimensional M = 1 
supersymmetric field theory is determined by the Kahler potential Kt(£ p , (pi) 
and the superpotential W(<pi), where £ p (p = 1, . . . , h 1 ' 1 ) represent the Kahler 
moduli and fa (i = 1, . . . h 2,1 + 1) the dilaton and complex structure moduli. The 
tree-level Kahler potential takes the form 

K T {£ p ,fa) = k{Q +K(fa), (5) 

where K is the Kahler potential for the Kahler moduli and K that for the dilaton 
and complex structure moduli. We will write an explicit expression for K in the 
next section. As for K, the only information we need to provide is that it is of 
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the no-scale form, in the sense that K~ l p q K p K q = 3 where p and q label the 
Kahler moduli. 

The appropriate superpotential for non-vanishing fluxes was proposed by 
Gukov, Vafa and Witten [22] and takes the form 

W = f (F 3 -T# 3 )Aa (6) 

J M 

t is the complex dilaton-axion field and Q the holomorphic (3,0) form for the 
Calabi-Yau M. 

The combination of the fact that W does not depend on the Kahler moduli 
and the no-scale structure of the Kahler potential implies that the effective scalar 
potential is simply given by 

V = e KT (d % WD 3 WK^j , (7) 

where DiW = diW + (diK)W (in Planck mass units). Therefore, to stabilise 
these fields at the minimum of the potential we need D^W = 0. The Kahler 
moduli only appear in the potential in the e Kr overall factor and are not fixed. 
We will omit these fields in the rest of this paper except for pointing out that the 
no-scale structure implies that the order parameter for supersymmetry breaking 
is 

D P W = (d p K)W. (8) 

Thus a non-vanishing value of the superpotential W implies supersymmetry 
breaking by the F-term of the Kahler structure moduli. Since a Kahler transfor- 
mation K — > K + / + / , W — > e~^W leaves the action invariant, the value of W 
can always be rescaled and a more appropriate, Kahler transformation invariant, 
measure of supersymmetry breaking is the gravitino mass, given by 

«3/2 = e K \W\ 2 . (9) 

Notice that in the KKLT construction, a non-perturbative superpotential de- 
pending on these fields is added |5|HT|. breaking the no-scale structure and fixing 
the Kahler moduli. Anti D3-branes (or magnetic fluxes on the D7-branes [T!Z] ) 
are introduced in order to lift the minimum of the potential to a positive value. 
This can certainly be done in the present model, but as we have nothing new to 
say in this respect we concentrate on the dependence of the theory on the dilaton 
and complex structure moduli, to which we now turn. 

2.2 Fluxes, Periods and Moduli 

For any Calabi-Yau 3-fold, the middle homology and cohomology are naturally 
expressed in terms of a symplectic basis. That is, there exists a basis of 3-cycles 
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A a and B b and a basis of 3-forms a a and f3 b (where a,b — 1,2... (/i 2,1 + 1)), such 
that in homology 

A a f]B b = -B b f]A a = 5%, 
A a n A b = B a n B b = 0, (10) 

and 



<*a = - P" = <t (11) 
A b JB a 



a a A{3 b = - (3 b Aa a = 5 b a . (12) 

M JM 

Such a symplectic basis is only defined up to Sp(2n, Z) transformations, as these 
preserve the symplectic intersection form. The periods are defined as the integral 
of the holomorphic 3-form Q over these cycles, 



n = z a } / n = g a . (13) 

A a JB b 

The periods are encapsulated in the period vector, II = (Qi, . . . , Q n , z%, . . . , z n ), 
where n = h 2,1 + 1. This inherits the holomorphic freedom of Q and is defined 
up to holomorphic rescalings Q — > f(<pi)Q. Note that II = 11(0$), where 0« are 
the complex structure moduli of the Calabi-Yau. 

Given the vector of periods II the Kahler potential on complex structure 
moduli space is given by |2S1 



K(r,<f>i) = -ln(-i(r-f)) -In I ij OAQ 

= -ln(-2(r-f))-ln(-2n t -S-n) (14) 



where 

V _ ( !n 

L "l-i„ o 

and r is the dilaton-axion. We can then compute the metric on moduli space, 

9ap = d a dpK, (15) 
and the Riemann and Ricci curvatures 

R^up — -dp(g Xa dp9pa), 

Rpp = RU = -d^\og(detg). (16) 
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In terms of the periods, the Gukov-Vafa-Witten superpotential is 



W=(2n) 2 a'(f-rh)-U( ( p i ), (17) 

where / = (fx, . . . , J 6 ) and h = (hi, ... , h 6 ) are integral vectors of fluxes along 
the cycles. The fluxes _F 3 and if 3 are elements of H 3 (A4,Z) and satisfy the flux 
quantisation conditions 



Fez, — — — — - / H ez 



(2ir) 2 a' y s3 ' (2Ti) 2 a' 



where S 3 is an arbitrary 3-cycle. The amount of D3-brane charge carried by the 
fluxes is then 

N flux = f T -Z-h. (19) 

The stabilisation of the complex structure moduli fa is governed by the fol- 
lowing equations 

D^W = d^W + {d^K)W = 0, 
D T W = d T W + (d T K)W = 0. (20) 

Distinct choices of fluxes stabilise the complex structure moduli at discrete points 
in moduli space. As the total number of flux vacua is very large, the distribution 
of discrete flux vacua can be approximated by a continuous distribution. Treating 
the fluxes as non-quantised, Ashok and Douglas ^E] derived a formula for an index 
that gives an estimate of the total number of vacua with Nf{ ux < L on a region 
X of moduli space 

WAW < L) = ^rj'^ 2 / det(-K -1-lo) (21) 

where u is the Kahler form and TZ the curvature two-form on moduli space. In 
formula (}2"Tj) K is the number of cycles along which flux is wrapped and L the 
total available D3-brane charge. In a comparison with numbers of explicit vacua, 
the constant prefactor is obviously not relevant. 

We will now restrict to cases where the Calabi-Yau Ai has two complex struc- 
ture moduli. The region X includes the dilaton-axion moduli space. The weighted 
vacuum density over the Calabi-Yau moduli space is then evaluated to be[28j 

dfi = g Tf dr A df A ( 4tc 2 c 2 - det^)^ 1 A dip 1 A dip 2 A dip 2 ) , (22) 



where g Tf = — 7—^2 an d c 2 , the second Chern class of M., is given by 

c 2 = — - (tr(ft Ml)- txTZ A trK) . (23) 
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m ay be rewritten as 



dfi = g Tf dr A df A c^ 1 A c^ 1 A dip 2 A eft/; 2 e°V b (R l a - al R 2 bl2 - R\- a2 R 2 bll ) - det # c 

(24) 

To evaluate both f!2()l) and (|22j) we must obtain a knowledge the periods, which 
in general is a highly non-trivial task. However, the manifold defined by (0) is 
of a class that has been extensively studied. The relevant periods have been 
computed in [201 El], following the classic treatment of the quintic[19 j; and we 
will borrow their results. For the Calabi-Yau described by the hypersurface 

4 d 

P = x d - k] — dif)x XiX2X 3 x 4: (f>xl° xl 1 xf xl 3 x\ 4 = 0, (25) 

3=0 yU 

the fundamental period in the large ijj region is given by 
where 

«,(#) - W)< £ ^ifa^Mz^!! . (27) 

S (' - D«)!"!nf =1 (^n + i(,„ - q,)l)\ 

This is obtained by direct integration of Q and satisfies the Picard-Fuchs equation. 
There are other independent solutions to the Picard-Fuchs equation having a 
logarithimic dependence on In total there are six independent solutions, one 
for each 3-cycle, and the actual periods are a linear combination of these. 

The two regions of moduli space that will most interest us are the vicinities of 
the Landau- Ginzburg point ip = (j) = and the the conifold locus 864?/) 6 + 0=1. 
To obtain a basis of periods in the small x/j region, we analytically continue (|2*Bj) 
to obtain 

d^r(n)r(i i))r(i- ^ )r(i-4f)r(i-^p) • 

(28) 

Here u u {4>) is related to the hypergeometric functions and is defined through the 
contour integral 

uM = ~ r^=(^-Cr- (29) 

The contour integral is initially defined for lm(0) > and then defined over the 
rest of the plane by deforming the integral contour. The branch cuts, which are 
unavoidable when v is non-integral, start at ±1 and run to ±oo. 
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We may derive a basis of periods from the fundamental period in a simple 
manner. If we define 



■(2ttz) 

~4T 



-w f (a j <ip,a jqo <j)) 



(30) 



then w(ip,(j)) = (<U7 (4>A), ™i(4>, 4>) , ^(ip , <f>) , m(4>, 0),^4(^,0), w 5 (tjj,(j))) gives 
a basis of periods known as the Picard-Fuchs basis. Naively (|30|) would seem to 
give 12 independent periods, but there are interrelations discussed in j2Hj- The 
net result is that, as expected, there are six independent periods. This basis is 
however not symplectic. A symplectic basis is given by 



n(-0, 0) = m ■ vc(ip, (j)). 
Here m is computed in [23 EHj and is given by 

/ 



(31) 
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In principle, equation ()31|) completely determines the periods near ip = 0. 
However, it involves the integral expression ()29|) for u u (<f)) which is inconvenient 
for a computational treatment. Such a treatment is facilitated by the power series 
expansion of u v {4>) in the small region found in [20]. This is given by 



u v (0) 



ITYV . 

e 2 



r 1 + 



2T (-u) 



E 

m=0 



e mm / 2 r (^p) (20)' 



m!r (1 - n^hx) 



(32) 



The region of convergence of equations (|2~Kj) . (j2Hjl . and (|3^jl is worth discussing. 
We will here be specific to the particular Calabi-Yau described in section |21 This 
manifold develops a conifold singularity when + 864'0 6 = ±1 and there is also 
a strong coupling singularity when = ±1. The regions of convergence are 
determined by the singularities. All three equations are only valid for | 1 < 1 
and equation (|3^|l has the additional restriction |0| < 1. 

We will also be interested in the periods near the conifold locus. As will 
be further discussed in section ()3.2|) . the periods here have a certain standard 
form. However, for their exact numerical determination, we will use the neolithic 
approach, directly evaluating the power series (J3~T|) near the conifold locus. 

Finally, the Calabi-Yau we work on is the original manifold M. defined by 
the locus of the polynomial (j2J), and not its mirror W. A4 has a total of 128 
complex structure moduli. We expect some to be removed by the orientifold 
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symmetry, but there will still be many which we are ignoring. The validity of 
this was explained in [21]. The group G is a symmetry of (J2J) and if we only turn 
on fluxes invariant under this symmetry then the superpotential can only have 
a higher-order dependence on the other moduli. It is thus consistent to set all 
other moduli equal to zero and focus only on the moduli in equation (j2J and their 
associated fluxes. We will comment briefly on the general situation in the last 
section. 



3 Vacua Statistics 

The two natural regions for testing the Ashok-Douglas formula are the vicinities 
of the Landau- Ginzburg point -0 = = and the conifold locus 864^ 6 + 0=1. 

3.1 The Vicinity of ^ = <j> = 

A symplectic basis for the periods was given in equation (J3T|) . Let us untangle 
this in the vicinity of -0 = 0. We can expand n(-0, 0) as 

n = a(0) + 6(0)^ 2 + c(0)^ 4 + O (^ 6 ) • (33) 

Here a, b and c are vector functions of whose explicit form arises from the 
combination of equations (J2EJ), (JH2J), (|3U|) and (|3*T|). It can be checked that a} ■ E ■ 
b = a} ■ S • c = 0, implying 

n 1 " • £ • n = (a f • S • a) + {b} ■ £ • 6)V^ 2 + C (|Vf ) , 

and consequently 

K^<f>) = -in(-int-E-n) 

= -ln(-mt-S-a) -In f 1 + p| ' S ' ~^ 2 + O 

V (o T • S ■ a) 

= - In (Hat • E • a) — ffi ' ^ *W + O (|^| 6 ) . (34) 
Equations (ffilj) then have the form 

1/^ = 0^ (/-r/i)-(« 1 (0) + a 2 (0)^ 2 + a 3 (0)^ 2 ) =0, (35) 

D^W = 0^ (f-rh)-^ 1 (cf>)+l 2 ((f>)^ = °> ( 36 ) 

D r V7 = (/ - fh) ■ (a(0) + b((P)ij 2 ) = 0, (37) 

where we have dropped terms of Cd^j 4 ). Here a(0) and /3(0) are complicated 
functions of depending on the integral expressions for u v {(j)). However, when 
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|0| < 1 the use of the power series expansion in equation coverts a(<f>) and 
(3{4>) to a more tractable form. The leading behaviour of the metric g a a = d a dgK 
is given by 

This is consistent with our expectations - at the Landau-Ginzburg point i[) = <fi = 
the metric becomes singular. The curvature 2-form 7£ and the Chern class C2 
may be calculated straightforwardly from the full expressions for the metric using 
a symbolic algebra program. Evaluating the Ashok-Douglas density, we find it 
has leading behaviour 

dfi ~ g Tf dr A df A ip^dip A dip A d(f> A cf0. 

The regions on which we compare the Ashok-Douglas formula to our empirical 
results are balls in and space. We then expect as leading behaviour 

iV(vacua s.t. < r x ) ~ r l5 
iV(vacua s.t. |0| < r 2 ) ~ rj. 

To test this expectation, we generated random fluxes and sought solutions 
of equations (J35|) to (pTTjl using a numerical root finder. The range of fluxes 
used was (-20, 20). This is not as large as one might prefer. However, a larger 
range of fluxes resulted in solutions being produced insufficiently rapidly for our 
purposes. In order to be able to trust our truncation of the power series, we only 
kept solutions satisfying 

< 0.5 and |0| < 0.75. (38) 

When processing the numerical results, there is an important subtlety we must 
account for 4 . It is well known that there is an exact SL(2, Z) symmetry of type 
IIB, 

ar + b ( F 3 \ f a b \ f F 3 \ 
T ^ cr + d 1 \H 3 J^\c d J\H 3 J' 

where a, b, c, d E Z and ad— be = 1. Thus each vacuum found has many physically 
equivalent SL(2, Z) copies that we should not double-count. One way to deal 
with this would be to fix the gauge explicitly and then perform the Monte-Carlo 
analysis. Our approach is instead to weight each vacuum by the inverse of the 
number of copies it has within the sampled flux range. The purpose of this is to 
ensure that vacua with many SL(2, Z) copies are not preferred. 

As well as the SX(2, Z) symmetry, there is a monodromy near the Landau- 
Ginzburg point that needs similar treatment. From the definition of the periods 

4 We thank S. Kachru for bringing this to our attention. 



864^ 6 
0± 1 
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(|3*0|). we can that they have a monodromy under (■0,0) — ► (a"0, — 0)> where 
a 12 = 1, of 

0) — > a ■ zu(ip, 0), 

where 

/ 1 \ 

1 

1 
a ~ 1 0' 

1 
\ -1 00000/ 

In the symplectic basis, the monodromy matrix A is given by m ■ a ■ mT 1 . The 
effect of this monodromy is to produce a family of physically equivalent solutions 
related by 

(0,0) -> (a~V,-0)> 

/ -> /-A 

-> fr-A (39) 

When weighting vacua we need to find the total number of copies lying within the 
sampled flux range from all symmetries and monodromies. This has important 
systematic effects as vacua with smaller values of fi and hi, and thus smaller 
values of Nfi ux , have more copies. A naive counting that neglects the symmetries 
or monodromies that are present thus places undue emphasis on vacua with 
smaller values of Nfi ux . 

We looked at the distribution of vacua within fixed balls in -0 and space. 
In figure ^ we plot the number of vacua satisfying (|3*K|) and having < r. The 
results are seen to agree well with the theoretical prediction. Likewise, figure El 
shows the distribution of vacua for a ball |0| < r in space. The continuous 
line again represents the cumulative number of vacua and the dots the rescaled 
numerical integration of J^ <r dfi. The empirical results are again well captured 
by the theoretical prediction. Finally, in figure El we examine the dependence of 
the number of vacua on the distance in flux space iVfl ux = f T ■ X • h. The graph is 
fit by iV ~ L 4 3 . This is surprising, as the expected scaling is L 6 . Furthermore, in 
the vicinity of the Landau- Ginzburg point for an analogous one-modulus example, 
the correct L 4 scaling is found 29 . As we will discuss further in section we 
believe our results are an artifact of the small flux range used, and that were a 
larger flux range used we would obtain the correct scaling. As we will shortly 
describe, in the vicinity of the conifold locus we do obtain the expected scaling 
with a flux range of (—40,40). 

It is also of interest to study the supersymmetry breaking scale, as measured 
by a' 2 (27i) A e K \W\ 2 , for vacua in the vicinity of = = 0. This is shown in 
figure |U The most noticeable thing about this graph is that the distribution of 
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1200 h 




Figure 1: Number of vacua with < r. The value of N plotted includes a 
weighting due to the SX(2,Z) copies of each vacuum within the range of fluxes 
sampled. The dots represent the numerical results, and the continuous line the 
(rescaled) numerical integration of J^ <r dfi, where dfi is the index density. The 

flux range used was (-20,20). 

the susy breaking scale is uniform near the origin and that the vast majority of 
vacua therefore have a high supersymmetry breaking scale. We can also observe 
from the structure of the graph that arbitrarily small values of the gravitino 
mass should exist as well as supersymmetric solutions. In reference [21], explicit 
super symmetric vacua were obtained for this model. 

3.2 The Vicinity of the Conifold Locus 

The Calabi-Yau M. has a codimension one conifold degeneration along the moduli 
space locus 864 , 6 + = 1. The moduli space curvature diverges near a conifold 
point and the expectation is that vacua should cluster near the conifold locus. 

To study these vacua numerically, we must restrict attention to a small region 
of the conifold locus where we can compute the periods explicitly. We take 
this region to be the neighbourhood of the point = 0, i\) = ipo = 864 _ e. If 
we write if) = i/; + £ and truncate the periods at first order in £ and 0, then 
II(£, 0) = (Qi, Qz, G3, Zi, z 2 , z 3 ) is approximated by 

Qx = 3202£ + 171.80 + C(£ 2 ,0 2 ,£0) 

G 2 = 4323 -z(1553£ + 107.40) + O(£ 2 ,0 2 ,£0) 

g 3 = (-492.7 + 1976.8z) + (371.0 - 300.2i)£ + (-259.0 - 59.Oz)0 + C(£ 2 , 2 , £0) 

Zx = ^C?iln((?i) + 784.8i-23O6z£-44.35z0 + O(£ 2 ,0 2 ,£0) 
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Figure 2: The weighted number of vacua, N, with |0| < r. The dots represent the 
numerical results and the continuous line the numerical integration of Jj^| <r d[i, 
rescaled by the same factor as in diagram Q 

z 2 = (-994.6 - 184.8?) + (861.9 + 476.5i)^+ (9.91- 112. 7z)0 + C(^ 2 ,0 2 ,^) 
z 3 = z(369.5 - 953.O^ + 225.40) + 0(£ 2 ,0 2 ,£0) (40) 

The numerical values were found by evaluating the series (|2*5Jl up to one hundred 
terms in ip and twenty-five terms in cj). The values for the coefficients of the 
0(£, <fi) terms were sensitive to the number of terms used in the power series at 
the level of a couple of percentage points. We did not keep terms quadratic in £ 
and cj) - inclusion of these would lead to 0(£) corrections to the results below. 

The general form of the periods is standard. The cycles corresponding to 
Q%t Gzi z 2 and z 3 are all remote from the conifold singularity, resulting in the 
associated periods being both regular and non-vanishing near the conifold degen- 
eration. Recalling that the conifold is a cone with a base that is topologically 
S 2 x 5" 3 , the cycle corresponding to G\ is the S 3 that shrinks to zero size along 
the conifold locus. The period along this cycle is regular near the conifold locus 
and vanishing along it. Finally, the cycle corresponding to Z\ is that dual to 
the S 3 . It is therefore not uniquely defined and is in fact well known to have a 
monodromy under a loop around the conifold locus. Its period takes the form 
(— 2^l^i hi({?i) + analytic terms). 

It is convenient to define Z = £ + <fi and rewrite the periods as 

Qx = 171.8 Z, 

g 2 = 4323 + 107.4zZ- 3554z£, 

G 3 = 784.8 - (259 + 59i) Z + (5198 + 799i) f , 
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Figure 3: The weighted number of vacua N with Nfi ux < L. 4.5 x 10 6 sets of 
flxues were generated, with values of L equally distributed between and 800 
and the range of fluxes being (—20, 20). The results are fit by iV ~ L 4 3 . 

Zl = ^171.8 Z In (Z) + 784.8 - UAi Z - 1479zf, 

z 2 = (-994.6 -184.8f) + (9.91- 112.7?) Z+ (677 + 2577z) 

z 3 = i(369.5 + 225.4 Z- 51540- (41) 

Z is a measure of the distance from the conifold locus. We are interested in 
vacua extremely close to the conifold locus - typically hi|Z| < —5 - and thus 
we will regard \Z\ « |£| << 1. Having set up the periods (}4~Tj) . we can now 
compute the Ashok-Douglas expectation for the index density and compare this 
with numerical results. 

Let us now solve equations (f271jl . First, 

D T W = (f - fh) • II = r = (42) 
This can be written as 

r = f±^l + 0(Z In Z\ (43) 
fro + &if 

where and &j are flux- dependent quantities. Next, 

D{W = (f - rh) ■ (c + d£ + c 2 l + O(e)) = 0. (44) 

Using (|4*3*|) this becomes a linear equation for £, easily solved to determine £ and 
r. We finally need the value of In Z. This is obtained by considering 

D Z W = => (f 4 - rh 4 ) \nZ= (d + d x r) + (d 2 + d 3 r)£ + (d 4 + d 5 r)£ (45) 
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Figure 4: The value of (27r) A e K \W\ 2 in units of (a') 2 for vacua in the vicinity of 
i/> = <p = 0, for {2n) A e K \W\ 2 < 100000. The flux range was (-20, 20), and the 
vacua satisfied the conditions (jHEJ)- 



Substituting in for r and £ from ()43|) and ()44j) then gives the value of In Z. 

When analysing the results we must account for the SL(2, Z) copies discussed 
in section IH.ll There is also a monodromy near the conifold. When solving for 
In Z, we impose no restriction on the imaginary part of In Z. There is then a 
monodromy 

\nZ -> \nZ + 2ni, 

(hi hi hi hi hi h) (h + hi hi hi hi hi h)i 
{h 1 ,h 2 ,h 3 ,h^h 5 ,h 6 ) -> (h t + h 4l h2,h 3 ,h 4 ,h 5 ,h 6 ). (46) 

corresponding to a loop in moduli space around the conifold locus. This gives 
a further source of physically equivalent solutions that should not be double- 
counted. For each vacuum, we then count the total number of copies within the 
specified flux range and weight by the inverse of this number. 

Results are shown below. In figure 03 we show the clustering of vacua by 
plotting the distribution of vacua transverse to the conifold locus. A similar plot 
of vacua parallel to the conifold locus shows no such clustering, indicating that 
there is no preferred position along the conifold locus. In figure El we perform 
a quantitative comparison with the expected vacuum density, finding very good 
agreement over a large range of values of In Z. In figure [7| we examine the scaling 
with L of the number of vacua having Nfi ux < L. This reproduces the expected 
L 6 scaling of (|2*T|) . We also computed the susy breaking scale, a' 2 (27v) 4 e K \W\ 2 . 
The distribution of susy breaking scales is plotted in figure El and is uniform near 
zero. 
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Figure 5: The value of Z for vacua near the conifold. We have restricted to 
\Z\ < 0.0001 and have rescaled the above plot by 10 4 . The flux range used was 
(-40,40). 

4 Summary and Discussion 

Let us briefly summarise our results. 

1. We constructed a large class of flux vacua for the two- moduli Calabi-Yau 
threefold : 2 2 6 j . We independently computed the Ashok-Douglas density 
and compared with our results. We find good agreement which improves as 
the range of fluxes is increased. The number of vacua was limited mostly 
by working in two patches in moduli space: the region near ip = <fi = and 
the region close to the conifold singularity <fi + 864-0 6 = ±1. 

2. We found a large concentration of vacua close to the conifold singularity 
as predicted, with the detailed distribution of the vacua being in close 
accordance with expectation. 

3. In both regions, the values of the superpotential are uniformly distributed 
near zero. As a consequence, large values of the superpotential, correspond- 
ing to a large supersymmetry breaking scale, are more abundant than small 
ones. 

4. In both regions the number of models scaled as a power of the maximum 
permitted value of Nfi ux , L max . In the vicinity of the conifold we reproduced 
the expected L 6 scaling. In the region close to the Landau- Ginzburg point 
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Figure 6: The distribution of vacua transverse to the conifold. We plot the 
number of vacua having InZ < D for D G (—120, —5) against D, restricting to 
|£| < 0.05. The dots refer to the vacua found numerically and the smooth line 
to the Ashok-Douglas prediction. We include results for three flux ranges - (-20, 
20), (-30, 30) and (-40, 40). The fit of the results to the expected distribution 
improves markedly as the flux range is increased. 

we only achieved an scaling as L 4,3 . We attribute the failure to achieve the 
expected scaling in this region to the smaller range of fluxes used there. 

An important feature of our analysis is that the flux range used has a sig- 
nificant effect on the results. If the flux range is insufficiently large, then the 
distribution of vacua found numerically will not fit with the theoretical density. 
This is most strikingly illustrated for the case of the conifold in figure |3 Similar 
behaviour was seen for the scaling of the number of vacua with L - in the vicinity 
of the conifold locus, a reduction of the flux range to (-20,20) caused a reduction 
in the power of the scaling from w 6 to ~ 5. Given this, we believe that the scal- 
ing of iV(vacua \Nji ux < L) ~ L 4,3 found in the vicinity of the Landau- Ginzburg 
point is simply an artefact of the flux range used. 

We can see a similar dependence on the range of fluxes used in the distribution 
of the dilaton. After being transformed to the SL(2, Z) fundamental region, the 
expectation from equation (|22jl is that the number of vacua with Im(r) > tq 
should scale like — . In figure 9 we compare this with numerical results arising 
from using flux ranges (-20,20), (-40,40) and (-60,60). 

We see that as the flux range increases the empirical distribution tends to- 
wards the theoretical one, and also that even with a flux range of (—60, 60) the 
two distributions do not yet fully match. In general terms these results are reas- 
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Figure 7: The weighted number of vacua with iVfi ux < L. The curve is fit well by 
N oc L 6 . The range of fluxes used was (-40,40). 



suring in the sense that arbitrarily large values of the dilaton can be obtained, 
consistent with the weak coupling aproximation although they are not statisti- 
cally preferred. 

As we have concentrated on only two complex structure moduli from a total of 
128, we have explored only a minuscule sample of the full set of flux vacua for this 
model. Still, we obtained sufficient statistics to see the expected distribution of 
vacua. There is a clear computational obstruction to exploring the full spectrum 
of models for a given Calabi-Yau. However, we have seen that the consistency 
requirements of small string coupling and large volume are largely disfavoured 
statistically. The relevant class of vacua will then be suppressed compared to the 
large estimated totals of order 10 100 . Nonetheless, given this large number it is 
reasonable to expect many solutions with very small values of the supersymmetry 
breaking scale (which in the KKLT scenario becomes the large volume constraint). 
It would be desirable to develop techniques that, rather than simply counting 
the number of vacua by numerically scanning the parameter space, select only 
those which satisfy the consistency and phenomenological requirements. Some 
small values of \W\ can be found explicitly - contains an example where 
\W\ «4x 10~ 3 and the first reference of jH2] a case where \W\ 2 « 2 x 10~ 4 . 
The efficient discovery of phenomenologically acceptable vacua is harder; the 
techniques recently developed in [3U] may be of use in this regard. Finally, let us 
note that in our study of this model we have addressed neither the issue of soft 
supersymmetry breaking terms induced solely by the fluxes nor a' corrections 
nor the non-perturbative superpotential required to complete the moduli fixing 
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Figure 8: The value of (27r) 4 e^|VF| 2 in units of (a') 2 for vacua near the conifold, 
for (2n) 4 e K \W\ 2 < 100000. The flux range was (-40, 40) and we restricted to 
vacua satisfying |£| < 0.05. 



and lift the potential to get de Sitter space as in the KKLT scenario. For recent 
progress in theses directions, [HU El may be consulted. 

Let us finish with a general observation. The Calabi-Yau considered here has 
h 1 ' 1 = 2 and h 2,1 = 128. We have turned on fluxes along the cycles corresponding 
to only two of the 128 complex structure moduli. We would expect that some 
of the remaining moduli should be frozen out by the orientifold symmetry, but 
it would still be obviously impractical to attempt either to write down or to 
solve the moduli stabilisation equations with all fluxes turned on. However, if we 
assume that the Ashok-Douglas density remains valid then we can say something 
about the generic situation. Suppose we have K cycles supporting flux and that 
- as holds for this and many other F-theory models - 

N D3 + N flux = L max ~ 1000. (47) 
The Ashok-Douglas density (|2~Tj) tells us that 

iV(vacua \N flux < L«) ~ Lf . (48) 

The fraction of vacua having Nd 3 > n is then estimated by 

iV(vacua \N fiux < L max - n) _ (L max - n) K 
N (vacua \N flux < L 

max ) ^max 

For K ~ 200, then for n = 5 this is approximately - « 0.36. We then see 
that despite the large amount of D3-brane charge we have to play with, generic 



21 



2000 



1000 



1500 



500 




4 



6 



8 



10 



Figure 9: Distribution of vacua near the conifold locus with Im(r) > y for y e 
(2, 10). Results have been brought to the same scale and plotted for three separate 
flux ranges: (-20,20), (-40,40) and (-60,60). We see that as the flux range increases 
the empirical plot moves closer to the expected result (represented by a smooth 
line) . 

vacua have rather small gauge groups 5 . This is appealing from a phenomeno- 
logical point of view if the standard model were to live on D3 branes as in the 
realistic models of [31]. Of course, the actual numbers depend on the particular 
Calabi-Yau, the number of 3-cycles surviving the orientifold projection and the 
modifications to the Ashok-Douglas density when L w (a few) if, but the general 
conclusion should be broadly unaffected. Clearly, there is no reason to require 
that our universe belong to statistically preferred classes of vacua, but given the 
discussion above, it is always useful to know which properties are generic and 
which are not. 
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